SLAC-PUB-13517 



Matter Dependence of the Three-Loop Soft Anomalous Dimension Matrix 

Lance J. Dixon 

SLAC National Accelerator Laboratory, Stanford University, Stanford, CA 94309, USA 

(Dated: January 27, 2009) 

The resummation of soft gluon exchange for QCD hard scattering requires a matrix of anomalous 
dimensions, which has been computed through two loops. The two-loop matrix is proportional 
to the one-loop matrix. Recently there have been proposals that this proportionality extends to 
higher loops. One can test such proposals by computing the dependence of this matrix on the 
matter content in a generic gauge theory. It is shown that for the matter-dependent part the 
proportionality extends to three loops for arbitrary massless processes. 
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High-energy colliders, particularly hadron colliders 
such as the Tevatron and the imminent Large Hadron 
Collider (LHC), copiously produce complex events con- 
taining multiple jets of hadrons. The production rate 
for such events can be understood through perturbative 
quantum chromodynamics (QCD). However, in many 
cases large logarithms of ratios of kinematic quantities 
appear, and it is necessary to resum these logarithms to 
obtain an accurate prediction. For processes in which 
only two partons appear at the leading order in the per- 
turbation expansion in a s (Born level), such as quark- 
antiquark annihilation to an electroweak vector boson, 
or gluon-gluon fusion into the Higgs boson, such resum- 
mations have been carried out to relatively high order [l[ , 
thanks in part to a detailed understanding of properties 
of the Sudakov form factor in QCD [1, y, H|. In these 
cases, only one color structure is permitted — the colors 
of the two partons must be identical. For three partons 
as well, a unique color structure appears at Born level 
(an SU(3) generator matrix or structure constant). 

However, for four or more partons, such as appear in 
di-jet production at hadron colliders, the Born level pro- 
cess contains multiple color structures. The emission or 
virtual exchange of soft gluons can mix these structures. 
To organize the logarithms from virtual exchange it is 
convenient to introduce a soft for cross ) anomalous di- 
mension matrix @, 0, [H, [l0| . This matrix can be 
computed from the renormalization properties of vacuum 
matrix elements of products of Wilson lines, also known 
as eikonal lines. Each line corresponds to an external 
parton in the amplitude. At one loop, the soft matrix 
can only change the colors of two hard partons at a time, 
because the virtual gluon propagator has only two ends. 
Such action can be written as a color operator of the form 
T°T° = T, ■ Tj [lUIli where T; indicates the action on 
the color of line i, and the adjoint index a of the virtual 
gluon is summed over. 

At higher loops, it might seem that arbitrarily many 
external partons (eikonal lines) could eventually be con- 
nected by soft gluons, resulting in a very complicated 
soft anomalous dimension matrix. For example, at two 
loops the structure f abc T[^ 
connecting eikonal lines i, j 
however, such a structure does not appear ^for any two- 
loop amplitude in a massless gauge theory [13|, [l4| ■ Also 



the two-loop soft anomalous dimension matrix T s • 
proportional to the one-loop matrix, 
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T^T£ might have appeared, 
and k. Rather remarkably, 
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where 7^ is the L-loop coefficient of the cusp anoma- 
lous dimension @], 1k{o. s ) = J2T=i Ik \ a sh) 1 



' ' — 2Ci (Cj is the quadratic Casimir for line i) 



with 
This 



result agrees with explicit infrared singularities of various 
two-loop QCD amplitudes [H, [l|| [l7f . The soft anoma- 
lous dimension matrix also appears in the study [l8j and 
resummation [T^ | of electroweak Sudakov logarithms, rel- 
evant for very high-transverse-momcntum processes at 
the LHC. Recently it has been conjectured that the pro- 
portionality <JI} might persist to higher loops 21, 22, 23]. 
Ref. [23[ also showed that such an ansatz is consistent 
with a required invariance under rescaling of Wilson-line 
velocities, but the solution to the consistency conditions 
is not unique for amplitudes with four or more partons. 
In this letter the proportionality (TI| is shown to ex- 

(3) 

tend to three loops, for the part of T s that depends 
on the matter content, in any massless gauge theory. 
The result is established using the same transforma- 
tion of loop-integration variables employed at two loops 
in refs. [H, [ijj]. There is also a term in Yg arising 
from pure gluon exchange. Although the result here 
shows that certain pure-glue contributions connecting 
three eikonal lines also vanish, other contributions con- 
necting three and four lines remain to be understood. 
On the other hand, these contributions are the same 
in any theory, so it is possible, for example, to extract 
them from an amplitude in J\f = 4 super- Yang-Mills the- 
ory (SYM), and then apply the result to QCD. In fact, 
the full color dependence of the three-loop four-point am- 
plitude in M = 4 SYM is known in terms of a handful of 
loop integrals [2l[. Evaluating the integrals (a nontriv- 
ial task) would allow a test of the proportionality of the 
pure-gauge part of the soft anomalous dimension matrix 
at three loops. The observed uniform trans cendental- 

ity of Af = 4 SYM amplitudes [24| already implies that 

(3) 

any non-proportional term in T s should be a weight-5 
transcendcntality function of kinematic invariants, where 
Q(n), In™ x and Li n (a;) are examples of weight n functions 
— in addition to the constraints of ref. [23| . 



Consider the scattering amplitude for n massless par- 
tons /, carrying (all-outgoing) momenta {pi} and color 
{r{}. To analyze the amplitude at fixed angles and large 
overall momentum scale Q, represent the momenta as 
Pi = Qv{, v 2 = 0, where the Vi are four- velocities. Di- 
vergences are regulated using dimensional regularization 
with D = 4 — 2e. The amplitude's color-dependence can 
be represented as a vector \M) in a C-dimensional vector 
space spanned by the color tensors {(c/)r r 1} 0, [H, HH , 
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For example, for n-gluon amplitudes, one choice of 
basis is the set of multiple traces of generators, 
tr(T Ql . . . T a ") tr(T afc +! . . . T a ') . . .. For amplitudes with 
two quarks (1 and 3) and two anti-quarks (2 and 4), the 
basis is much simpler; it is spanned by S^S^ and S^S 1 ^. 

Infrared divergences of on-shell amplitudes may be 
factorized into jet functions describing virtual partons 
collinear with the external lines, and soft functions char- 
acterizing the exchange of soft gluons between the hard 
lines. The general form of the factorized amplitude, for 
equal factorization and renormalization scales /i, is [251 ] 

n 

|X(^,Q 2 //i 2 ,a s ( M ),e)) = JJjM(a a 0*),e) 

i=l 

xS(v,,Q 2 / t i 2 ,a s ( f i),e)\H(v ]7 Q 2 / f i 2 ,a s ( t i))) , (3) 

where is the jet function for external state i, S is the 
soft function, and H is the hard (short-distance) func- 
tion, which is finite in the infrared. 

The jet function for parton i can be expressed to all 
orders in terms of three anomalous dimensions, K\ l \ 
and 7^; and /C^ is completely determined by 7^ 0, |3]- 
This letter concerns the soft function S. It is governed 
by the soft anomalous dimension matrix Tg, 



S (vi ■ Vj Q 2 //i 2 ,a s (/i),e) 
dfl T 
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A* 2 
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and also involves the /^-dimensional running coupling. 
At one loop, a s is given by 
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with p = (11C A - 4T fl n/)/3 in QCD. 
The one-loop soft matrix is 0, 



4 X) = JeEt s .t, 



In 
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FIG. 1: Three- loop diagrams containing a fermion or scalar 
loop and three eikonal lines, corresponding to a three-gluon 
vertex at two loops. 
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FIG. 2: Three-loop three-eikonal diagram containing a 
fermion or scalar bubble, and not originating from a three- 
gluon vertex at two loops. 



where Sij = (pi +pj) 2 . Resummed cross sections are 
determined by the eigenvalues and eigenvectors of T$ [t| 



[10|,l26|,l27|. Proportionality of T^ 2) toT^ implies that for 
fixed kinematics, the same eigenvectors that diagonalize 
the one-loop matrix also diagonalize the two-loop matrix, 
simplifying the integration over scale [l3|, [l4[ ■ 

The soft anomalous dimension matrix is found from 
the single-ultraviolet-pole terms in e, for suitable combi- 
nations of eikonal lines, as described in detail at one Q 
and two loops [lil ]. At L loops, the maximum number 
of eikonal lines that can be connected to some gluon is 
2L. However, this maximal number corresponds to a dis- 
connected product of L one-loop single gluon exchanges, 
which appears in the expansion of the exponential ([4]). 
Thus the divergence is removed by one-loop renormal- 
ization. The maximum number of eikonal lines that does 
not correspond to a disconnected product of lower-loop 
configurations is L+l. Requiring dependence on the mat- 
ter content of the theory means that there is a fermion or 
scalar loop in the diagram; this cuts the maximum num- 
ber of nontrivially connected legs to L. Such diagrams 
can be generated from the maximal nontrivial (L — 1)- 
loop configurations by inserting bubbles with fermions or 
scalars, or triangles, or higher-point one-loop diagrams. 

At two loops, the key eikonal diagrams for demonstrat- 
ing proportionality involved three eikonal lines [l3l . [l4| , 



,(i) 



3 



labelled by A, B and C, with velocities v A , vb and vc- 
The diagrams can be divided into two cases: (i) three 
eikonal lines connected uniquely by a three-gluon vertex, 
and (ii) three eikonal lines connected by two single-gluon 
exchanges, so that two gluons attach to one of the eikonal 
lines. Diagram (i) was found to vanish. The sum of di- 
agrams of type (ii) was found to have only a symmetric 
color configuration, proportional to {T^ • Tc , Tc -Tb}; 
the coefficient factorized into an appropriate form to be 
removed by one-loop renormalization. (It corresponds 
to second-order terms in the expansion of the exponen- 
tial (UJ) in which one leg is common to both factors.) 

Inserting a matter triangle or bubble into diagram (i) 
gives the three-loop diagrams in fig.QJa) and (b). Insert- 
ing a bubble into a type (ii) diagram gives diagrams of the 



form shown in fig. [21 I shall show that this insertion does 
not affect the properties established for the pure-gluon 
diagrams at two loops. The proof uses exactly the same 
change of integration variables used at two loops, plus, 
for diagram [TJa), a handy representation of the one-loop 
three-gluon vertex / aia2a3 r /JlM2M3 (fa, k 2 , k 3 ), evaluated 
in an arbitrary gauge [28l |. Inserting this representation 
into diagram fjja) yields the loop-momentum integral, 

F 1(a) (v A ,v B ,v c ) (x f abc T A T b B T c 

j 1D1 dPk \_VAVBVG /Q\ 

J 1 2 k\k\k\v A - h v B -k 2 vc ■ k 3 ' 
omitting the +ie prescription, with k 3 — —fa — k 2 and 



^v A v B v c = "A u fl"c r mM2»( fc i. fc 2, fa) = A(k\, k\; k\)v A ■ v B v c • (fa - fa) + B(k\, k\\ k\) v A -v B v c ■ (fa + k 2 ) 

- C(kl,kl;kl)(k 1 ■ k 2 v A ■ v B - v A ■ k 2 v B ■ fa) v c ■ (fa - k 2 ) 
+ ^S(kl, fcf, fcf ) (v A ■ k 2 v B ■ fa v c ■ ki + v A ■ fa v B ■ fa vc ■ k 2 ) 

+ F(k\, k\; kl) (fa ■ k 2 v A ■ v B - v A ■ k 2 v B ■ fa) (k 2 ■ k 3 v c • fa - fa • k 3 v c ■ k 2 ) 

- H(k\,kl,kl)[v A ■ v B (k 2 ■ k 3 v c ■ fa - fa ■ k 3 v c ■ k 2 ) — ±(v A • k 2 v B ■ k 3 v c ■ fa - v A ■ k 3 v B ■ fa Vc ■ k 2 )] 
+ cyclic permutations of (fa, A), (k 2l B), (k 3 ,C). (7) 

I 



The precise forms of the coefficient functions A through 
H are given in ref. 28], for gluon and massless fermion 
loops. Here only their symmetry properties are needed: 
A, C and F are symmetric under exchange of the first 
two arguments; B is antisymmetric in the first two argu- 
ments; S is totally antisymmetric; and H is totally sym- 
metric. Ref. [28[ does not explicitly consider the case of a 
massless scalar, but the vanishing of three-point functions 
in M — 4 SYM implies that the scalar loop must have the 
same tensor decomposition and symmetry properties. 

The change of loop-momentum variables used in 
refs. [H],[lII was expressed in terms of light-cone variables 
using the vectors v A and vb to define + and — directions. 
Here it suffices to give the action of the transformation 
on the various Lorentz products that appear, 

V A ' Vc j . , 

va - fa —> vb • ki, vc ■ ki — > vc ■ ki , 

V B ■ V C 

vb ■ h — > VB VC v a ■ fa , fcj • kj — > fa ■ fcs , (8) 
v A ■ v c 

where 1 = 2, 2 = 1, and 3 = 3. The Jacobian for this 
transformation is unity. It is simple to check that ev- 
ery contracted tensor appearing in eq. (|7|) has the proper 
symmetry under |8]) so that, combined with the symme- 
try of the appropriate coefficient function, and the re- 
maining (symmetric) product of propagators in eq. 
the integrand is antisymmetric. Hence the integral ^ 



vanishes. 1 For example, the two terms multiplying S in 
eq. ([7]) exchange with each other with a positive sign, 
whereas the corresponding terms multiplying H do so 
with a negative sign. Note that the transformation (JSJ) is 
only to be used for the term shown explicitly in eq. ijTj). 
For the two images of this term under cyclic permutation, 
one should use the correspondingly permuted transfor- 
mation, which respects the symmetry properties of the 
permuted functions A through H. 

For the diagrams in fig. QJb) an d fig- [21 write the bubble 
diagram as 

n£i%( fc2 ) = * aiQ2 (w 2 - ^A 2 A 2 )n(fc 2 ) , (9) 

where £ = 1 for the transverse scalar and fermion bub- 
bles, but £ may differ from 1, depending on the gauge, for 
a gluon bubble. The ?7 MlAt2 term in eq. ©, inserted into 
diagrams |Tfb) and [21 clearly gives the same tensor struc- 
ture as at two loops. Thus [13|] its contribution to dia- 
gram[ljb) vanishes. It is easy to see that the ^-dependent 
terms are also odd under the transformation ©, so dia- 
gram [IJb) vanishes for any £. 

There are four diagrams of the type shown in fig. [21 re- 
lated by changing the order of connection on line C , and 



Another argument for the vanishing uses the total antisymmetry 
of the integrated result under exchanges of {va , Vb , vc }, plus 
invariance under suitable velocity rescalings 12911 , 



4 



by moving the matter bubble to the other gluon line. 
The sum of the four diagrams splits into a color antisym- 
metric piece, proportional to [T^ • Tc , Tc ■ Tb], whose 
coefficient vanishes exactly as in ref. [l4[; and a color 
symmetric piece, proportional to {T^ • Tc , Tc • Tg}, 
containing a factor in the integrand of 



1 



va ■ vc I v B ■ v c - £ 
k 2 



vc ■ h 

+ (va <-» V B , h 
v c ■ (fa + k 2 ) 



, , v A ■ v c [v B ■ v c - £ 

vc ■ k x v c ■ k 2 [ \ 

x n(fcj) + (v A ^v B , k 1 ^k 2 ) 



+ (va <-> v B , k\ f-> k 2 ) 
vb ■ hv c ■ ki 



kl 



(10) 



The latter form implies that this diagram factorizes into 
the product of a one-loop diagram with the matter- 
dependent part of the two-loop diagram. Hence it also 

(3) 

does not contribute to T y s . The above arguments also 
show that, independent of the gauge, the soft matrix does 
not receive contributions from those three-eikonal pure- 
glue diagrams containing a gluon loop. 

Finally one must address the matter-dependent con- 

(3) 

tributions to T s involving only two eikonal lines. As 

(2) 

was the case for the full T y s ' , only two of these dia- 
grams (ladder and crossed ladder) do not obviously re- 
duce group theoretically to the form T^ • Tj . Inserting a 
matter bubble of the form ([9]) does not change the group 
theory. The magnitudes of these two contributions are 
linked [2(| l30| . however, by the general non-abelian ex- 
ponentiation theorem for eikonal webs [3lj , in such a way 



that their sum is of the form T^ ■ T^ (after removing the 
square of the one- loop term). Because this result holds 
for the case of two partons (i and j of equal colors, so 
that T.; • Tj — > T^ = Ci), the constant of proportion- 
ality must be the (matter dependent part of) the cusp 
anomalous dimension, 



i(3) , matter 



Tk 



77 



A' 



,(1) 

■ s 



(11) 



An alternative, very general argument for the constant 
of proportionality, given that the matrix is propor- 
tional, is based on an anomaly under rescalings of the 
eikonal- line velocities (23|. As noted there, the situa- 
tion may become more complex at four loops, due to 
the fact that 7^ mi ght contain quartic Casimir terms 
- (df cd ) 2 [la [M 111. ( At ver y stron g coupling, 
"quadratic Casimir scaling" is known to break down for 
planar Af = 4 SYM (H.) 

In summary, the results of this letter show that the 
three-loop soft anomalous dimension matrix is essentially 

the same in all massless gauge theories, i.e. up to the 

(3) 

matter dependence of "f K ■ If proportionality fails for 
the pure gauge terms, then at least the three-loop prop- 
erties are constrained by JV = 4 super- Yang-Mills theory, 
as well as by the velocity rescalings studied in ref. [23] . 
Acknowledgments 

This work was supported by the Department of Energy 
under contract DE-AC02-76SF00515. I thank Adi Ar- 
moni, Thomas Becher, Zvi Bern, Gregory Korchemsky, 
George Sterman, and especially Einan Gardi and Lorenzo 
Magnea, for helpful correspondence. 



[1] S. Moch, J. A. M. Vermaser en and A. Vogt, Nucl. 
Phys. B 726, 317 (2005) liep-ph/0506288] , and refer- 
ences therein. 

[2] A. H. Mueller, Phys. Rev. D 20, 2037 (1979); 
J. C. Collins, Phys. Rev. D 22, 1478 (1980); J. C. Collins 
and D. E. Soper, Nucl. Phys. B 193, 381 (1981) [Err.- 
ibid. B 213, 545 (1983)]; A. Sen, Phys. Rev. D 24, 3281 
(1981). 

[3] J. C. Collins , Adv. Ser. D irect. High Energy Phys. 5, 573 

(1989) [hep-ph/0312"336] , 

[4] L. Magnea and G. Sterman, Phys. Rev. D 42, 4222 

(1990) . 

[5] R. A. Brandt, A. Gocksch, M. A. Sato and F. Neri, Phys. 

Rev. D 26, 3611 (1982). 
[6] A. Sen, Phys. Rev. D 28, 860 (1983); J. Botts and 

G. Sterman, Nucl. Phys. B 325, 62 (1989). 
[7] G. P. Korchemsky and A. V. Radyushkin, Phys. Lett. 

B 171 (1986) 459; S. V. Ivanov, G. P. Korchemsky and 

A. V. Radyushkin, Yad. Fiz. 44 (1986) 230 [Sov. J. Nucl. 

Phys. 44 (1986) 145]. 
[8] N. Kidonakis, G. Oderda and G. S terman, Nucl. Phys. B 

531, 365 (1998) [hep-ph/9803"24i] . 
[9] N. Kidonakis, G. Oderda and G. Sterman, Nucl. Phys. B 



525, 299 (1998) [hep-ph/9801268] . 
[10] A. Banfi, G. P. Sala m and G . Zanderighi, JHEP 0408, 

062 (2004) [hep-ph/0407287] . 
[11] S. Cata ni and M. H. Se ymour, Phys. Lett. B 378, 287 

(1996) [hep-ph/9602"277| ; Nucl. Phys . B 485, 29 1 (1997) 

[Err.-ibid. B 510, 503 (1998)] [hep-ph/9605"323] . 
[12] S. Catani, Ph ys. Lett. B 427, 161 (1998) 

[hep-ph/9802439] , 
[13] S. M. Aybat, L. J. Dixon and G. Sterma n, Phys. Rev. 

Lett. 97, 072001 (2006) [hep-ph/0606254] . 
[14] S. M. Aybat, L. J. Dixon and G. Ste rman, Phys. Rev. D 

74, 074004 (2006) (hep-ph/0607309| . 
[15] Z. Ber n A. De Freitas an d L. J. Dixon, JHEP 0306, 028 

(2003) [hep-ph/0304"i"68] . 
[16] E. W~ N. G lover, JHEP 0404, 021 (2004) 

[hep-ph/0401119] . 
[17] Z. Bern, L J. Dixon and D. A. Kosower, JHEP 0408, 

012 (2004) [hep-ph/04042931 . 
[18] B. Jantzen, J. H. Kiihn, A. A. Penin and V. A. Smirnov, 

Nucl. P hys. B 731, 188 (2005) [Err.-ibid. B 752, 327 

(2006)] [hep-ph/0509T57l . 
[19] J.-y. Chiu, R. Kelley and A. V. Manohar, Phys. Rev. D 

78, 073006 (2008) [0806.1240 [hep-ph]]. 



5 



[20] I. A. Korchemskaya and G. P. Korchemsky, Nucl. Phys. 

B 437, 127 (1995) [hep-ph/9409446] . 
[21] Z. Bern, J. J. M. Carrasco, L. J. Dixon, H. Johansson and 

R. Roiban, Phys. Rev. D 78, 105019 (2008) [0808.4112 

[hep-th]]. 

[22] T. Becher and M. Neubert, 0901.0722 [hep-ph]. 

[23] E. Gardi and L. Magnea, 0901.1091 [hep-ph]. 

[24] A. V. Kotikov and L. N. Lipatov, Nucl. Phys. B 661, 19 
(2003) [Err.-ibid. B 685, 405 (2004)] [hep-ph/0208220] ; 
Z. Bern, M. Czakon, L. J. Dixon, D. ~K. Kosower 
and V. A. Smirnov, Phys. Rev. D 75, 085010 (2007) 
[hep-th/061 0248 ; S. G. Naculich, H. Nastase and 
H. J. Schnitzer, JHEP 0811, 018 (2008) [0809.0376 [hep- 
th]]. 

[25] G. Sterman an d M. E. T ejeda-Yeomans, Phys. Lett. B 

552, 48 (2003) [hep-ph/0210"l30] . 
[26] R. Bonciani, S. Catani, M. L. Mangano and P. Nason, 

Phys. Lett. B 575, 268 (2003) hep-ph/0307035 . 



[27] Yu. L. Dokshitzer and G. Marchesini, Phys. Lett. B 631, 
118 (2005) [hep-ph/0508l"30] ; JHEP 0601, 007 (2006) 
[hep-ph /0509078 j; M. H . Seymour, JHEP 0510, 029 
(2005) hep-ph/0508305 ; A. Kyrieleis and M. H. Sey- 
mour, JHEP 0601, 085 (2006) [hep-ph/05 10089] , 

[28] A. I. Davydychev, P. Osland and O. V. Tarasov, Phys. 
Rev. D 54, 4087 (1996) [Err.-ibid. D 59, 109901 (1999)] 
[hep-ph/9605348j . 

[29] E. Gardi, private communication. 

[30] G. P. Korchemsky and A. V. Radyushkin, Nucl. Phys. B 

283, 342 (1987). 
[31] J. G. M. Gatheral, Phys. Lett. B 133, 90 (1983); 

J. Frenkel and J. C. Taylor, Nucl. Phys. B 246, 231 

(1984). 

[32] L. F. Alday and J. M. Maldacena, JHEP 0711, 019 

(2007) [0708.0672 [hep-th]]. 
[33] A. Armoni, JHEP 0611, 009 (2006) [hep-th/060 8026 . 



